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Abstract 



A new weak measurement procedure is introduced for finite samples which 
£T) . yields accurate weak values that are outside the range of eigenvalues and which 

do not require an exponentially rare ensemble. This procedure provides a 
£^ ■ unique advantage in the amplification of small non-random signals by mini- 

[ — . mizing uncertainties in determining the weak value and by minimizing sample 

size. This procedure can also extend the strength of the coupling between the 
^ system and measuring device to a new regime. 

Ck 
-i— > ■ 

1 Introduction 

Aharonov, Bergmann and Lebowitz (ABL, pQ) considered measurement situations 
between two successive ideal measurements where the transition from a pre-selected 
^ ■ state l^in) to a post-selected state |^fi n ) is generally disturbed by an intermediate 

precise measurement. A subsequent theoretical development arising out of the ABL 
work was the introduction of the "Weak Value" (WV) of an observable which was 
probed by a new type of quantum measurement called the "Weak Measurement" 
(WM) [3] (reviewed in £J2]). The motivation behind these measurements was to explore 
the relationship between |^i n ) and |\l/fi n ) by reducing the disturbance on the system 
during the intermediate time. This can be modeled by reducing the interaction 
strength between the system and the measuring device. For example, if a WM of 
A is performed at the intermediate time t (t- m < t < tg n ) then, in contrast to the 
ABL situation, the basic object in the entire interval t- m — > for the purpose of 
calculating other WVs for other measurements is the pair of states I^V) an d l^fin)- 
However, the reduction of disturbance also reduced the information obtained from 
a single WM on a single quantum system. Therefore, the WV was determined by 
using a large ensemble (reviewed in §2.11 and §2.21) . This was a result of the weakness 
condition which produced a shift in the pointer of the measuring device (MD) that 
was much less than the uncertainty. Many separate irreversible recordings of the 
slight MD shift were then used to amplify the "weak value signal" above the "noise" 
due to the weakened measurement. 
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This article introduces a new Gedanken experiment coined "Robust Weak 
Measurements on Finite Samples" (RWM). The primary advantage of RSM is a 
reduction in the uncertainty of the WV for finite samples and also increases the 
probability to obtain WVs which are outside the eigenvalue spectrum. RWM involves 
an irreversible recording of the sum of momenta for an ensemble of quantum systems, 
such that the shift in the sum of momenta is large compared to it's noise. In addition, 
in order to ascertain the maximally allowed information about the WV, the relative 
positions (which commute with the total momenta) are also measured. In a physical, 
"realistic," WM, there is always a finite coupling and thus a disturbance caused to 
the system in addition to unknown fluctuations. However, we can use the relative 
positions to correct for this disturbance and unknown fluctuations and thus WVs can 
now be determined much more accurately for a finite-sized ensemble. We illustrate 
RWM by a practical application of WVs to the amplification of weak or unknown 
signals and show how the interaction strength A can be increased and yet still have 
a useful regime of WVs. RWM also allows for a reduction in the number of particles 
necessary to perform an accurate WM. 

2 Weak Measurements 

WMs [31 [11] can be quantified in the quantum measurement theory developed by 
von Neumann [IT] . First we consider an ideal measurement of observable A by using 
an interaction Hamiltonian Hi nt of the form H mt = — X5(t)Q md A where Q md is an 
observable of the MD (e.g. the position of the pointer), A is a coupling constant which 
determines the strength of the measurement, and 5(t) determines the duration of the 
measurement (setting % = 1). For an impulsive measurement we need the coupling 
to be strong and the duration short and thus take 8(t) to be non-zero only for a short 
time around the moment of interest such that J T 8(t)dt = 1 (which thereby allows 
us to ignore the free Hamiltonians of the system and MD). Using the Heisenberg 
equations of motion for the momentum P md of MD (conjugate to the position Q m d)> 
we see that P md changes according to = X5(t)A. Integrating this, we see that 
Pmd(T)— P m d(0) = XA, where P m d(0) characterizes the initial state of MD and P m d(P) 
characterizes the final. To make a more precise determination of A requires that the 
shift in P m d, i.e. <5P m d = P m d(P) — P m d(0), be distinguishable from it's uncertainty, 
AP m d. This occurs, e.g., if P m d(0) and P m d(T) are more precisely defined and/or 
if A is sufficiently large. However, under these conditions (e.g. if MD approaches a 
delta function in P m d), the disturbance or back-reaction on the system is increased 
due to a larger H mt which is a result of the larger AQ m d arising from the inverse 
relationship of AP m d and AQ md (AQmd > irk — )• When A is measured in this way, 

then any operator O ([A, O] ^ 0) is disturbed because 40 = iX5(t)[A, 0}Q md , and 
since AAQmd is not zero, O changes in an uncertain way proportional to AAQmd- 
In the Schroedinger picture of measurement, the system and MD state are: 

\$m) = \^in)K d ) -> e~ 4 / M |v|/ in >|<^ d > = e a4 ^ in )fc> (1) 

where the state of the system is |^i n ) and the MD state, l^md); i s given by |$™ d ) = 
f dQ md $™ d (Q md )\Q md ) = J dP<&^ D (P m d) | P m d) • A good approximation for realistic 
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experiments is to consider MD's initial state as a Gaussian (without loss of gen- 

erahty), e.g. $£ d (Q md ) = (Q md fc> = exp(-gf) and $^(P md ) = (Pfc) = 
exp(— A 2 P 2 ) (substituting A = AQ m d, A? m d = i, leaving off the normalizations). 
Expanding \^m) in eigenstates of A, i.e. \^m) = J2% \A = di)(A = ai\^/ in ) = 
J2i a i\A = Oj) then eq. [U becomes: 

X></ rfQmde a ^e-&|A = a,>|Q md > = E^/ dPe^ A ^ A V = a,)|P) (2) 



When the uncertainty, AP md = , in MD is much smaller than the shift of MD 



5Pmd — ^ a h corresponding to the strength of interaction and to different eigenvalues, 
then the final state of MD is a density matrix representing a series of peaks, each 
corresponding to a different eigenvalue a*, i.e. after tracing over the state of the mea- 
sured system, the absolute square of eq. [2] yields Pr(P md ) = J2i Wi\ 2 e~ 2A i p ^-^i) _ 
In other words, MD goes into a state of superposition proportional to the system. 
If 3>^(-P m d ~ ^ a i) is orthogonal to $j^(P md — Xaj) when % ^ j (which occurs if 
<5P m d AP md ), then this macroscopic superposition collapses into a single peak 
with probability given by the Born rule Pr(aj,t\^i n ,ti n ) = \ (aj[Ut in ->t\^m) | 2 = | Q j| 2 , 
i.e. depending only on the initial state of the measured system. 

However, suppose at a later time tg n , we measure another non-degenerate operator 
which yields, as one of its outcomes, the state l^fin). The conditional probability to 
obtain a,-, given both boundary conditions, the pre-selected \^m) and post-selected 
(tffij, is given by ABL 0: 



Pr(aj,t\^ in ,t in ; ^ &n ,U n ) 



Prj^fin, hn\aj, t)Pr(aj, t|^ in , t L 
E tt ' Pr(*fin,tfin|a / ,t; V- m ,t) 







^in)| 2 






H 


*in)| 2 



(3) 



To illustrate the time-symmetry in ABL, we apply U t ^t &n on (^fij, i-e. ( l I / fi n |Pt_ > t fln = 
(C/ t Ut fin *fm| and use UUt^ = {e""^-*)}* = e iH ^~^ = e -«*(*-fc») = [7^^. The 
time-symmetric re- formulation of the numerator of eq. EJ (Ut &n ^t^B.n\o>j)(0'j\Ut in -*t\^hi), 
can now be interpreted to mean that the time displacement operator Ut Rn -*t sends 
(\IZfinl back in time to t as depicted in fig. [Ha. The Born rule, of course, is recovered 
by summing over all possible final states {^fin}^ 

Pr( aj , t\y in ; t in ) = J2 Pr(aj, t\V in , t in ; {* fin } n , t fin ) Pr({* fin } n , t &n \ aj , t; * in , t) . (4) 



We can weaken the interaction H int = — \5(t)Q m dA by minimizing AAQ md . For 
simplicity, we consider A 1 (assuming without lack of generality that the state of 
the MD is a Gaussian with spreads AP md = AQ m d = 1)- We may then set e~ lX ® mdA ~ 
1 — iXQmdA and use a theorem |19^l : 

A\V) = (A)\V)+AA\y ± ), (5) 

1 where (A) = |*) is any vector in Hilbert space, AA 2 = (V\(A - (A}) 2 \^), and 

is a state such that (^l^j.) = 0. Note that (A) is not defined here in a statistical sense: it is a 
mathematical property of an individual system l^S). 
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b) System Measuring Device 

£fin 



'Pmd' 

£P md = Aai>AP md 
Ideal "strong" measurement 



fin 




l^ir) 5P md = AA w <AP md 



Weak measurement 



Figure 1: a) with an ideal or "strong" measurement at t (characterized e.g. by <5P md = 
}a\ S> AP md ), then ABL gives the probability to obtain a collapse onto eigenstate a\ by 
propagating (fl/finl backwards in time from tfj n to t and j^in) forwards in time from t- m to 
t; in addition, the collapse caused by the ideal measurement at t creates a new boundary 
condition |ai)(ai| at time t, t Ui < t < ifi n ; b) if a WM is performed at t (characterized e.g. 
by <5-Pmd = AA^ <C AP md ), then the outcome of the WM, the WV, can be calculated by 
propagating the state (^finl backwards in time from ig n to t and the state fd/j^ forwards in 
time from t m to t; the WM does not cause a collapse and thus no new boundary condition 
is created at time t. 



to show that before the post-selection, the system state is: 

e' A4 " di |*in) = (l-^4ndi)|*i„) = (l-^gmd(i))|*in)-AgmdAi|* in± ) (6) 
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Using the norm of this state || (1 - i\Q m dA)\^ in } || = 1 + A Q md (A 2 ), the prob- 
ability to leave l^m) un-changed after the measurement is: 



1 + A 2 QL(^) 2 



1 (A-0) (7) 



while the probability to disturb the state (i.e. to obtain \^i n ±)) is: 



A 2 Q md A^ 

+ A 2 g md W 



(A 0) 



The probability for a collapse decreases as -h. Thus, for a sufficiently weak inter- 
action (e.g. A < 1), the probability for a collapse can be made arbitrarily small. 
In addition, the measurement becomes less precise because the shift in MD is much 
smaller than its uncertainty (i.e. SP md = Aaj <C AP md ). The final state of MD is now 
a superposition of many substantially overlapping Gaussians with the same distribu- 
tion as before Pr(P md ) = J2i \(A = Oil^m) | 2 exp j — ^^"'^ j. However, this ends 

up being a single Gaussian ^ d (P m d) « (P m d\e- iX ^ d{A) K d ) » exp j- {Pmd ~^ ))2 } 
centered on A (A). 

If we perform this measurement on a single particle, then, of course, we will not 
be able to distinguish between two states which are not orthogonal, e.g. |$j° d ) 
and |\l/ 2 )|$™ d ). Such an ability would violate unitarity because these states could 
time evolve into orthogonal states |\l>i)|$™ d ) — > |$™ d (l)) and |^ 2 )|*md) — > 
|* 2 )|$j n d (2)}, with |*i)|$Sd(l)> orthogonal to |* 2 > l*md( 2 ))- Fr o m a WM per- 
spective, the reason that this does not happen is that measurement of these two 
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non-orthogonal states causes a smaller shift in MD than it's uncertainty and there- 
fore we might conclude that the shift 5P m d of MD is a measurement error be- 
cause iC(^md) = (P md - A(A)|$£ d ) « (P md |$£ d ) for A < 1. Nevertheless, 
if a large (N > 4-) ensemble of particles is used, then the shift of all the MDs 
(SP^ ~ A(A) ; y- = N'(A)) can accumulate to a distinguishable level while the col- 
lapse probability still goes to zero. That is, for a large ensemble of particles which are 
all either |\I/ 2 ) or this measurement can distinguish between them even if \^> 2 ) 
and are not orthogonal because the scalar product (\E f ^ N ' ) |\l/2 N ' ) ) = cos n 6 — > 0. 
Traditionally, it was believed that if a measurement interaction is limited so there 
is no disturbance on the system, then no information is gained. However, we have 
shown that when considered as a limiting process, the disturbance goes to zero more 
quickly than the shift in MD and thus with a large enough ensemble, information 
can be obtained even though not even a single particle was disturbed. 

Now that we have a new measurement which does not cause a collapse, we ask 
whether this type of measurement might reveal new values. We shall see that with 
a WM (which involves adding a post-selection to this weakened von Neumann mea- 
surement), the MD registers a new value, the WV. As an indication of this, we insert 
a complete set of states {|^fi n )j} into (A): 



(A) = <sg 



If we interpret the states |^/fi n )j as the possible outcomes of a final ideal measurement 
on the system (i.e. a post-selection) and we perform a WM (e.g. with \AQ md — > 0) 
during the intermediate time t, t- m < t < t&a, then the coefficients | (^fiiilj^in) j 2 give 
the probabilities Pr(j) for obtaining a pre-selection of (^hj and a post-selection of 
|^ / fin)j (since the intermediate WM does not disturb these states) and the quantity 
A w (j) = ^~[^y is the WV of A given a particular final post-selection (\l/ fin \j. 

Thus, from (A) = Y^j Pf(J) A W (J), one can think of (A) for the whole ensemble as 
being built out of pre- and post-selected states in which the WV is multiplied by a 
probability for post-selection. 

To see how the WV arises from a weakened measurement with post-selection more 
precisely, we consider the final state of MD in the position representation: 



C D (Qmd) = (Q m d|(*fin|$ to i) = (* fi n|e a ^ di |vl> in )e-^ 

= e { ~ tXQ r d) \ ^\% n )e- Q -^ = <* fil j* in > e HA ^ md) V ) 

n=0 nl n=0 nl 

= (* fin |*in){l + iXQ^A* + E { ^7 dr K}e~^ 

— 71 ' 

n=2 ,L - 

= (* & r\%n){l + i\Q md A w + g - ^^^^ (A w )- 

— 71 ' 

n=2 n - 

n=2 n - n=2 H - 



The second term in the last part of eq. [TO]) can be neglected in 2 general regimes: 

1. minimizing AAQ m d by either using a small A (setting AP m d = AQmd = 1), or by 
minimizing the spreads in MD, e.g. so that P m d is measured to a finite precision 
AP m d, (which limits the disturbance by a finite amount AQ m d > or 

2. minimizing [A™ — (A w ) n ]/n! even if AAQ m( j is not small. 

By way of example, the first moment in the Taylor's expansion (from the second term 
in the last part of eq. [JU]) can be neglected if (AAQ m d) 2 AA w <C 1 where AA W = 
\(A 2 ) W — (A w ) 2 \2 [3]. When eq. [10] is transformed back to the P m d representation, 
then the final state of MD after WM and post-selection is (up to normalization): 

SCOPmd) « (^ fin |^ n )(Pje- iA4 - A -|^ d )^exp{-A 2 (P md -AA w ) 2 } (11) 

^ ere = (12) 

The final state of MD is almost un entangled with the system and has the same initial 
shape but shifted by a very surprising amount, the WV, A w (the factor (^fml^m) 
arises as a result of the exclusion of other post-selections) . Since the value of A is 
given by P md (T) — P m d(0), we may conclude A « A w . We have used such limited 
disturbance measurements to explore many paradoxes (see, e.g. [5]). There have also 
been a number of experiments to test the predictions made by the WM and their 
results are in very good agreement with theoretical predictions [UJ [9] [13], [18] EE]. 

The new developments presented in this article are motivated by an application 
of WMs and WVs to quantum metrology which can provide a unique advantage 
over the usual (pre-selected-only) approach in the amplification of small non-random 
signals. In the pre-selected-only approach, the outcome of quantum measurements 
are restricted to the eigenvalue spectrum range. The advantage of WVs is that they 
can be far outside this range due to the overlap (^finl^in) in the denominator of 
A w , eq. [T2] These WVs are termed "eccentric weak values" (EWV). Now, if the 
pre- and post-selection are known to high precision, then in the idealized weak limit 
(e.g. AAQ m d — > 0), the WV can be calculated precisely and one might conclude 
that no new information is obtained if a WM is actually performed since MD will 
simply register a shift by XA W . However, if the coupling A between system and 
MD is unknown and contains additional small errors which are not random, then 
actually performing a WM which yields an EWV can provide new information, e.g. 
by allowing us to distinguish the shift from the non-random force (incorporated into 
A) from the large EWV shift due to the WM interaction. 

The intention of this article is to address several theoretical issues concerning 
such an application. In particular, we note that the WV approximations presented 
above become more and more precise in the idealized weak limit of AAQmd — > and 
N — > oo in which there is no disturbance or back-reaction on the system. Neither 
of these limits are realistic in practical amplifications because first of all, we must 
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have a finite N, and second of all with a finite N, we must increase A to discern the 
WV from the noise, and therefore we can no longer ignore the two uncertainties in 
determining the WV which arise due to: 

1. the inability of MD to measure definite WVs due to the MD's uncertainty AP m( j 

2. the back-reaction on the system due to AQ m d creates an uncertainty in the 
WV of the system itself 

In this article, we demonstrate a new approach to WMs, coined "Robust Weak 
Measurements on Finite Samples" (RWM) which decrease these uncertainties and 
allow us to increase the coupling (e.g. A ~ 1, not A < 1 as used in eq. [11]) for 
finite N while maintaining the accuracy of the WV. In order to accomplish this, we 
will have to consider the structure of MD in Q md (in addition to P md used above). 
Intuitively, we can see 2 inverse roles for MD observable Q md and system observable 
A. We have already reviewed how Q m d generates translations in P md proportional 
to A w . However, the roles for Qmd and A are reversed when one asks the reverse 
question: 'What is the back reaction on the system (i.e. not on the MD) due to the 
measurement interaction H mt 7" For this question, A is now the generator (not Q m d) 
in a manner proportional to Q m d ( n ot A). The structure of MD in Q m d therefore 
registers information concerning the back reaction on the system. This is measured 
with RWM, thereby reducing the second uncertainty (which results in PPS-mixtures). 
The new RWM introduced in §3] uses components of both WM criterion. The first, 
minimizing XAQ m d, is introduced in greater depth in the next section and the second, 
minimizing AA W is introduced in greater depth in §2.21 



2.1 Statistical Weak Measurements (SWM); A > -7=— 

V TV A w 

To make RWM more "Gedanken-practical," we will consider Stern-Gerlach (SG) 
measurements in different contexts throughout this article where A will be a spin 
component (i.e. <7g) and Q md will be the translational coordinate of the particle in 
the same direction, £ as the spin component. Having the particles themselves serve 
as MDs with the information about the measurement stored in a degree of freedom 
(P md ) separate from the pre- or post-selection (so that there is no coupling between 
the variable in which the result of the measurement is stored and the post-selection 
device) provides the easiest way to pick out only those MDs which are associated 
with those systems that satisfied the proper post-selection criteria: the post-selection 
of the particles then also selects out the relevant MDs. 

Suppose we pre-select a spin-1/2 system with \a x = +1) = \*\ x ) at time t- m . To 
do this, we filter a x — — 1 out of the beam by applying an inhomogeneous magnetic 
field, described by H int = —fia x B x where B x = xB . The force on the particle 
= fiB a x induces a change in momentum proportional to the time T that the 
particle spends in the field, i.e. <5P md = fiB a x T. Since the particle is constrained 
to be in a region Ax < D (with D the size of the Stern-Gerlach opening), the 
initial uncertainty in the momentum must be AP md > For this pre-selection 
measurement to create a distinguishable split between a x = +1 and a x = —1, the 
shift induced in the momentum by the inhomogeneous magnetic field must be greater 
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than the uncertainty in the momentum, i.e. 5P^ d > AP^ d . We then perform a 
similar procedure to post-select \fr y = +1) = | f y ) at time £g n . 

If we now consider SG measurements in the intermediate time t, t- in < t < tg n 
at an angle £ to the x — y plane, then the spin fr^ can be determined from the 
deflection of the particle which is proportional to the impulse SP^ d = Xfr^ imparted 
to the particle due to the inhomogeneous magnetic field which has a linear gradient 
in the same direction £ as the spin component which is to be determined. Since 
the particle is free, the spin is conserved in time and thus a measurement of either 
fr x or by at t will yield +1. This is also evident from ABL: the probability to 
obtain erg = +1 in the intermediate time if an ideal measurement is performed is 
Pr(a ? = +1) = i+cos(0+^(|)+««(0 sin(g) _ Wg gee th&t if £ = Q o ^ e & j then the 

intermediate measurement will yield fr x = +1 with certainty and when £ = 90° (i.e. 
fry), then the intermediate measurement will again yield fr y = +1 with certainty. 
Consider measuring the spin in a direction £ = 45°: 

gg = fr x cos 45° + fry sin 45° = — ~t5 y (13) 

v2 

From the results Pr(fr x = +1) = 1 and Pr(fr y = +1) = 1, we might wonder if we 
could simply plug in both their values fr x = +1 and fr y = +1 into eq. [13] and obtain 
<7£ = = 73 = \/2- Such a result would obviously be incorrect for an ideal mea- 
surement because the eigenvalues of any spin operator, including fr^, are ±1. We can 
also see from ( Ux ^ v ) = a ^ +(T y + ' T ^' 7v+av ' Jx = = 1 but implementing the above 

argument, we expect ( ^jp ) ~ jlTf) = ^ 7^ 1- Performing this step of replacing 
fr x = +1 and a y = +1 in eq. [13] can only be done if fr x and fr y commute, which 
would allow both values to be simultaneously definite. The probability statements 
are only simultaneously true if we do not perform fr y before fr x , since this would de- 
stroy |^in) = |ta)- So, in general, the finding that fr x = +1 with certainty or fr y = +1 
with certainty in the pre- and post-selected ensemble only held when one of these 
two measurements was performed in the intermediate time, not both. The physical 
reason that a measurement of fr% doesn't produce v^2 is that the particle is exposed to 
a magnetic field with a strong gradient in the £ = 45° direction, which causes the spin 
to revolve around this axis in an uncertain fashion. In other words, the conditions 
for an ideal measurement SP^ d = Xfr^ ^> AP^ d will also necessitate AQ^d ^ ^ 
which will thereby create a back-reaction causing a precession in the spin such that 
A© ^> 1 (i.e. more than one revolution), thereby destroying the information that in 
the past we had fr x = +1, and in the future we will have fry = +1. 

However, there is a sense in which both Pr(fr x = +1) = 1 and Pr(fr y = +1) = 
1 are simultaneously relevant for measurements in the intermediate time: if the 
measurement of fr x and fr y is performed (i.e. when fr^ is measured) in such a way 
that measurement of one does not disturb the other, which is precisely what occurs 
in a WM. For such a WM, the inhomogeneity in the magnetic field induces a shift 
in momentum which is less than the uncertainty SP^ d < AP^ d and thus a wave 



2 The filter only interacts with the component of the spin which is not transmitted, e.g. the 
a x = — 1 component would receive a strong repulsive interaction (via a potential a x B x — B x using 
a homogenous B x ), while the a x = +1 component would not have any change in it's momentum. 
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packet corresponding to <7x ^ y = 1 will be broadly overlapping with the wave packet 

corresponding to ax ^ v = — 1 because the deflection 5Q\ ad oc (a^ ■ VB) will not be 

discernable from the noise, i.e. SQ^ d <C AQ^ d (where AQ^d i s the dispersion in 
the particle beam). Thus, it cannot be determined whether any individual particle 
corresponds to CTx ^ y = ±1- Furthermore, since the MD is quite imprecise, we 
cannot say whether the distribution of results in the pointer was due to the original 
uncertainty in the pointer XAQ md , or due to the distribution of the observable of the 
system being measured, AA. 



cither 



or 



-1) 



WM 

Of <7 ?= 4 5 o 



On, 



■1) tfin- 



"1)/ 



+ 1 



1) 



1) 



all \a x = +1} particle 1 ^ particle 2y particle 3 

Figure 2: Statistical WM ensemble. 



particle N 



Nevertheless, with statistical weak measurements (SWM), the WV can be ob- 
tained robustly in a statistical sense from the mean reading of many separate point- 
ers (see fig. [2]). For example, an ensemble of N separate systems and N separate 
MDs are used (which again for SG are independent degrees-of-freedom of the same 
particle being measured). For each individual system, in between it's pre- and post- 
selection, one of the N MDs weakly measures the observable A of this single system 
and the outcome of this measurement is individually recorded. This is repeated for 
each of the N different systems, each with a different MD. After the post-selection, 
the subset of those MDs which were associated with those systems which satisfied 
the post-selection criterion are collected out of the larger set of all possible post- 
selections (if all post-selections were included, then the decomposition, eq. [9] would 
be reproduced). While the WM was performed during the intermediate time t, we 
will obtain the same outcome and a simplified analysis if the readout of the WM 
MD is delayed until after the post-selection. A statistical analysis is then performed 
on the results of only those MDs associated with the proper post-selection and an 
average is manually calculated. This statistical procedure reduces the uncertainty in 
the mean position by -7=, thereby allowing for a more precise calculation of A w [3]. 
When thus correlated with the post-selection, the measurement result (which was 
confirmed experimentally for an analogous observable, the polarization [T3j) is: 



[0"£ = 4 5 O 



(T„ \^$t\U) _ {(T y \a y } + {a x \U}} _ (t» |1 + 1| U) 



y/2 (14) 



For an individual spin, the component of spin a? is an eigenvalue, ±1, but the WV 
(&t)w — \/2 is a/2 times bigger, (i.e. lies outside the range of eigenvalues of a ■ n) 
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and is thus called an "eccentric weak value" In this case, we can use EWVs 

to amplify either the magnetic moment fi or the gradient in the field. 

If we now consider a finite sample of particles, we ask what are the minimum 
number of particles necessary in order to distinguish the WV "signal" (<x§) w from 
the noise, i.e. such that the total momentum shift is greater than the total deviation 
fiPmd ^ A-^md . E.g. suppose the momentum is deposited onto a photographic plate 
after the post-selection. From the WM interaction with each of the N particles, 
the total momentum deposited on the plate will be 5P^J = NX(a^) w . Now, for 
simplicity, we set AP^ d = 1 and thus the dispersion is also AP^ d = (AP^ d ) 2 = 
1. The total dispersion is AP^ = ^AP^ d = NAP^ d . However, the standard 

deviation AP^ = y AP^ and therefore, in order to distinguish the outcomes, 

SPjM = N5Pi d = NX(a^ w > AP^ = VNAPl d = VN. Therefore, the center 
of this distribution can be determined by using an ensemble of iV particles, where 

What, generally, are the limitations of this method, SWM? Suppose N = 20 and 
that in order to satisfy 5Pj^ d ^ AP^ we are content with 5P£2 > 2AP I J ld ) . In this 
case, we must increase the coupling A to overcome the fluctuations, so from 5Pj^ d > 
2AP£? and with (<x ?=45 °)w = V2, we have X(a^) w N > 2^N or A > j^^f = 

J±. However, with the SWM requirement of minimizing XAQ md , a WM interaction 

strength of A = is too large for a valid WM. When such a WM is attempted, 
AAQmd will have a significant back-reaction and will thereby create uncertainty in 
the boundary conditions. The reason we are interested in the minimum of A rather 
than its maximum is that we need to obtain an EWV (i.e. outside the eigenvalue 
spectrum), rather than an ordinary WV, in order to implement our amplification 
scheme. Now that we are dealing with finite samples, there will of necessity be a 
back-reaction on the system due to the WM. If the back-reaction rotates the pre- or 
post-selection by too much (e.g. by tt) then we will not obtain an EWV but rather an 
ordinary WV. In addition, this will rotate the pre-selected state by an uncertain angle 
and therefore we will not know the WV with certainty as was required to implement 
the amplification scheme because we don't know whether the outcome we obtained 
is related more to /x or to A w . In other words, up to this point we have argued that 
as a limiting process (i.e. when there is no back- reaction), we can measure the WV 
that would have been there even if no measurement was actually made. However, 
the back-reaction invalidates this approach. So, to control the rotation of the pre- 
or post-selection we need to control XAQ which is responsible for the back-reaction. 
Thus, while we are certainly able to measure A w up to 0(^=), we do not know which 
A w we are actually measuring. 

This limitation of SWM (i.e. N > y. , or alternatively A > - ) is 

altered by RWM which reduces the uncertainty in A w that is created when the back- 
reaction of MD changes the post-selection in an uncertain manner. 



3 WVs even further outside the eigenvalue spectrum which therefore offer the possibility of even 
greater amplification can be obtained by post-selecting states which are more anti-parallel to the 
pre-selection: e.g. if we post-select the +1 eigenstate of (cosa)u x + (sina)cr 2 , then (<r z ) w = Atan 
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2.2 Single Trial Weak Measurements (STWM) 

The limitation of SWMs, A > > , is significant when dealing with finite sample 
sizes. We can do better, however, by measuring collective observables In addition, 
this allows us to measure all WVs with great precision in one single (though pre- 
viously thought to be rare) experiment. We consider again our theorem: A\*ff) = 
{A)\*f>) + AA|\I/_|_) [IS]. We can also measure this property with no reference to 
statistics by applying this identity to a composite, iV-particle state (which can also 
be viewed as a single system such as a large spin) |\I/( N )) = |^) 1 |^) 2 ....|^)at and using 
a "collective operator," A^ = Y!£Li Ai (where A\ is the same operator A acting 
on the z-th particle). Using this, we are able to obtain information on (A) without 
causing a collapse and thus without using a statistical approach because any product 
state |\I/( N )) becomes an eigenstate of the operator A^ N \ To see this, consider [HIE] 

A(N)|^(N)). 

iv(i)l^ (N) ) + AA^ri^L N) (z)) 



^(N)l^(Nh = 1_ 

1 ' N 



(15) 



where (A) is the average for any one particle and the states (i)) are mutually 
orthogonal and are given by |^ N) (z)) = |W) 1 |tf) 2 ...|tf L ) i ...|W) JV . That is, the ith state 
has particle % changed to an orthogonal state and all the other particles remain in 
the same state. If we further define a normalized state |^/j N ' ) ) = J2i "^vl^j N ' ) (*)) then 

the last term of eq. Hi is ^=|#i N) ) and it's size is |^=|#i N) )| 2 oc ± -> 0. Therefore, 

|\I/( N )) becomes an eigenstate of A^ N \ with value A, as N — > oo (the second term 
decreases as 0(N~ 1 ^ 2 ) even if the particles are not all in the same state, as long as 
the composite iV-particle state is a product state). 

We shall now consider a similar setup as used in §2.11 We perform a WM of the 
collective observable in the 45° angle to the x — y plane of er^ = J2iLi &\ |2j- Using 
H^i = -^QmdE^i^, a particular pre-selection of |t a ) (i.e. \^ ) = Iljli \Vx)j) 

and post-selection \] y ) (i.e. (^ } | = Uk=i(h\k = Un=i{\U)n + i\iz)n}), we will 
show that the pointer is robustly shifted by the the same WV obtained in §2.11 i.e. 
y/2: 

( ^ )w = V5tf(( T ,|T,»" = ^ ± <W (16) 



Using [crf^, o^] = j^eijkHn&ki we wm see that with N — > oo, all these opera- 
tors commute: 

i<^run>i = i<if^E3?>i = i^ek)i < 

we see that for any given state and sufficiently large N, we may neglect the fact 
that these operators do not commute. In addition, for sufficiently large N, we may 
measure &^ (or 6"( N )) and the probability for a collapse can be made arbitrarily 

small. Using again the theorem, we have ^l^i?) = (4 N) )I^£°> + A 4 N) l*in!)- 
We will now show that as N — > oo, — ► 0. In addition, from this we conclude 
that at the limit, l 1 ^^) is an eigenstate of which means that not even one of the 
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spins will collapse. To calculate Aaf) we need (<rf)) = (^En=i<^> = fEL(^> 
and because the spin states are identical (^ N ') = We will also need: 



<(4 N, ) 2 > 



now since <r" and a™ operate in different spaces, 

<(^ N) ) 2 > = ^ (iV((^) 2 ) + /V(/V - f)(^,) 2 ) = (alf + 1(((^) 2 ) - (alf) 

and finally, Aaf ) 2 = <(af >) 2 > - (a^f = 1«(^) 2 > - (alf) = 0(1) 

(17) 

To obtain eq. [HI let us now calculate the final state of MD after post-selection: 

N ( \ N \ N 

i*r> = n<Tvii«p iQ m dE4' n ir.)iio (is) 

3=1 l iV k=l J i=l 

Since the spins do not interact with each other, we can calculate one of the products 
and take the result to the Nth power and eq. [19] can be re- written: 

IC D > = fl(T,l,exp|lg md ^| \uh\^T) = {(T,|exp|lg md a 5 | |t,)}I^ D > 

(19) 

Using the following identity exp {iaa^} = cos a + ia^sina [2U|, this becomes: 



AQmd .ys . AQ 



cos — — ta^ sin 



md 



N s AT 



iV 



JV 



= [(UU] ] cos ^ ia w sm—^-\ |$ in ) (20) 

where we have substituted a w = (o^)™ = % ^ . . We consider only the second part 

(the first bracket, a number, can be neglected since it does not depend on Q and 
thus can only affect the normalization): 

(J> MD = J 1 — ^ 2 (^ md ) 2 _ ^""iQmd } |<T) MD \ (01) 

fin ~~ | i\r 2 JV J 11 

As /V -> oo, we use (1 + § )* = (1 + f )^ a « e a and obtain |$^ n D ) « e iAaTO ^w |$^ D ). 
When projected onto P m d, this results in the same shift \/2 as the SWM example in 
§2.11 The maximum of A can be increased up to e\fN from the weakness condition 
and eq. [TO the weak uncertainty A(<t^) w ~ In this case, an individual spin 

is acted on by U = exp j^Omd^f n which rotates the spin around the £-axis by an 
uncertain angle A6 pa AA ^ md = -j=. Thus, the probability that an individual spin is 
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still in its original state is 1 — ^ and the probability that the entire N spin system 
stays in it's original state is {1 — jf} N ~ exp{— e 2 } — ► 1 for e <C 1. Thus, as N — > oo 
the intermediate WM will give (<r^) w = ^T^n^ry^ = v^2 robustly since the shift in 
MD 5P m d = N\(a^)w is greater than the uncertainty in MD, AP md = 1. A single 
experiment is now sufficient to determine the WV with great precision and there is no 
longer any need to average over results obtained in multiple experiments as we did in 
the previous section. Therefore, if we repeat the experiment with different MDs, then 
each MD will show the very same WVs, up to an insignificant spread of -j= (assuming 
we obtain the particular, rare, post-selection). Therefore, the information from both 
boundary conditions, i.e. \^- m ) = Ilili |Tx)« an d (^fi n | = n£i(tj/|i> describes the 
entire interval of time between pre- and post-selection (for plots see [T5l [T6] ). 

However, as pointed out in the beginning of this section, STWM has a major 
advantage over SWM in the minimum required values for A. There is no difference 
in the momentum shifts for SWM and STWM as both cases are SPj^J = N(a^) w X. 
However, for STWM, AP£? ~ 1, whereas for SWM, AP^ ~ V^V. To obtain 
fiPmd ^ A-^md > the minimum value for the coupling constant A for STWM can be 
smaller (A > N ^ ) by a factor of ^= than for SWM (A > . ). There is thus 

a regime of A for which the amplification scheme is invalid for individual particles, 
but for which it is valid for collective observables therefore substantially increasing 
the utility of STWM. 

However, STWM has a major short-coming that is resolved by RWM. While 
STWM is a valuable Gedanken experiment, the probability K^fml^m)! 2 ^ for all N 
particles to end up in the same final state |^fi n ) becomes exponentially small. With 
the particular STWM considered in this section, we have N particles pre-selected 
with a x — 1, a WM of 045° (which doesn't significantly disturb the spins which are 
thus still in the state a x — 1 after the (T450 measurement) and followed by a post- 
selection in the y-direction. The probability to obtain a y = 1 is 1/2 and thus the 
total probability of finding all N spins with a y = 1 is an exponentially small 2~ N . 

On the other hand, SWM requires a much smaller sample. For the particular pre- 
and post-selection used in §2.1[ approximately N/2 out of N pre-selected particles 
will satisfy the post-selection criterion and thus this result is not a rare outcome, 
making it much more attractive for our amplification scheme. 



3 Robust Weak Measurements on Finite Samples 
(RWM) 

The RWM introduced in this section shares positive attributes of both the "statis- 
tical" ( §2.ip and "single-trial" ( §2.2p approaches and seeks to minimize the 2 uncer- 
tainties discussed in §2] which resulted from a finite sample size. 

1. We keep a significant benefit from the first SWM approach, namely the PPS- 
ensemble necessary to obtain EWVs was not rare. However, the SWM has a 
major disadvantage for finite samples, namely in order to distinguish the WV 
from the noise, a relatively large coupling constant is required, A > ^ 
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which can cause a significant back-reaction on the system and therefore increase 
the second uncertainty in A w , creating a mixture of WVs. 

2. In the SWM approach, we measured each individual shift 5P m d and obtained 
the same total shift 8P£2 = N(a^) w X as was obtained for STWM, which did 
not involve a measurement of "individual" shifts. Therefore, with SWM we 
are in a sense doing too much because there is additional information which 
can be obtained (namely the relative positions which commute with the total 
momentum) if we wait to make the measurement of SP^J as in STWM. We 
can then use the relative positions to correct for the disturbance caused to 
the system which resulted from a stronger measurement interaction which was 
required by the use of part of the SWM approach. 

Without loss of generality, we present the new RWM in the framework of SG 
measurements used in previous sections: consider a large collection of particles where 
the MD is simply the position and momentum of those particles (see fig. ED- We 
perform the following a) filter out | j x .) at time £ in ; b) perform a WM of ax ^ y at the 
intermediate time t but wait until after performing the post-selection to read out the 
result of the sum of ~ | of these interaction^; c) filter out {[ y \ at time tfi n ; d) absorb 

the particles onto a photographic plate and measure the sum of momenta J2i=i Pmd 
(without measuring the individual P^)] this recording will produce a definite shift 
by a WV; e) measure the relative positions to determine what the pre-selected and 
post-selected system the WM in step (b) was a measurement of. 



Figure 3: RWM by Stern-Gerlach apparatus for weakly measuring 045°. 
Let us consider again a particular pre-selection for each particle of If^) and post- 



selection in the state |t y ), so |* in ) = nj=i |T*)j and l*fln> = IIj=i |T»)j- The WM 



interaction in step b) [2j is described by an interaction Hamiltonian which couples 



<7 5 with Q md of MD i.e. H int = XS^Q^&l (where a\ = {^} 3 ). This generates 



4 While we present this Gedankenexperiment in the same spirit as the Einstein Gedankenexper- 
imcnt, we also recommend that a WM of a x uses a field B (a x x — a y y) with a small width in y. 
There is then little variation in the wavefunction in the y-dircction. The x-dircction would not be 
constrained and the wavefunction can freely vary in x. With this method, only a shift in the x 
direction would occur; the wavefunction in y would always be left in the ground state because the 
force is too small to excite it . 




'45 



e) measure relative positions 



14 



shifts in the individual momenta due to the WM interaction (as occurred in §2.11) . 
However, unlike §2.11 there is only one irreversible recording of the sum of these 
shifts (as occurred in §2.21) . i.e. one irreversible recording of the total momentum 
of the N particles which were deposited onto a single photographic plate, followed 
by N measurements of the positions (used to deduce the N — 1 relative positions). 
This state of the photographic plate after its interaction with the N particles shall 
be referred to as the "final state of the MD". After the WM interaction and post- 
selection (but before the irreversible recording), the final state of MD is: 



N 



l^rnd) = U(Us^*QL%}\U)s\*&i) (22) 
j=l 

Here we have set the coupling to each spin to be / XS(t)dt = A and without loss of 
generality have taken the initial state of MD as simply a Gaussian in the coordinate 

f (Q* ) 2 1 

Q md of each particle, i.e. ($™ d ) 3 = exp < — . (A ™f )2 \- As will be seen later, it will 

prove useful to reformulate the MD observables in terms of two complementary, non- 
commuting, collective observables, a pointer corresponding to the sum of momenta, 
Pjndi an d it's conjugate Q^l which generates shifts in the pointer P m N J: 

N pj 

rirS = E^f (23) 



j=l 

n ni 



These definitions are particularly useful because if the uncertainty in the individ- 
ual Qmd' s i s ^Qmd ~ 1> then the uncertainty in Q^l * s a ^ so ^Qml ~ 1 due to 

[Qmd, p md\ = 1 (The spread in E-Ii OL is VN and thus ^ j ^ md wl). Using 



_ (Qmd) 2 



■A ^Qld) 2 



4(AQL) 2 



N 

E {Qmd 
j=l 



n VJ N n (N) 



md 



[Q 



(N)i2 
mdJ 



N 



4(AQ J 



Q 



md 



N' 

(N) 

md |2 



j=l 



N 



N 



N' 



[Q 



(N)i2 
md J 



(25) 



'md) j=l 

it will also be useful to re-write the wavefunction of MD as: 



II exp 



(Qmd) 



4(AQ 



exp {- 



(QL N i) 2 , r UQL 



md) 



4 ( A Qmd) 2 



}exp{- 



4(AQ< N) ^ 2 



■} (26) 



< md , 



We now show that measuring the relative positions provides corrections to the 
pre- or post-selection, thus giving a different WV for each particle, represented by 



-,(N) 



V .{Q 3 ■ md ) 2 

a 3 m (which thereby explains the utility of exp { - — md j^y^- — })• 
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3.1 Use of relative positions 

Besides the sum of momenta, we can also measure the N — 1 relative positions 
(without disturbing the system), 

*i = <&d-£^f = Q ^-77W (27) 



This is because [Q md — E-jvSE-Pmd] = [^'-^md] = which is easy to see be- 
cause each pair of relative positions commutes with the sum of momenta, i.e. [Qmd — 

Qmd' E n =l Pmdi = [Qmd' E n =l -^md] — [Qmd' E n =l -^md] = ^~ * = ^' usm g [Qmd' E n =l -^md] 

[Qmd' ^md] = i- Furthermore Q md - Eli % = Eli 

As a preparation to obtain both a measurement of the relative positions and of 
the total momenta we re-write eq. [22] as: 

N N An N An 

= n(Tj^xp{iA{Q^-E^Ka ex p{ iA E^aitx.) J 

j=l n=l JV n=l JV 



(Qmd)L\ r E 3 (Q m --^ 2 



x exp{- -^d^ }exp{ 



)(N), 
md 

'md ^/jy 



4(AQ^) 2 4(AQ 



md ; 



AT g(N) 

IKtj/lj exp{iA%o-|} exp{iA— tx)j 



j=i 



(«) 2 e,(qL-# 2 



x exp{ ^ m *A }exp{ cnIn > ( 28 ) 

How is this re-formulation of eq. [22] in terms of the relative positions Xj = 
Qmd ~ E helpful? To see this, we'll consider eq. [28] one particle at a time. 
For the jth particle, we can apply the first exponential of eq. EH exp{iXxjd^}, 
to either the pre-selected state \1 x )j 01 to the post-selected state (since the 2 
exponentials commute). What does this exponential do to the pre- or post-selection? 
As mentioned in £j2] Q^l anc ^ ^ (or in this case Xj and at) have 2 inverse roles: 
the back reaction on the system is generated by <j| in a manner proportional to ctj. 
However, the be measured exactly and can thus be replaced by a number. 

Therefore, exp{iXxjd^} simply rotates the pre- or post-selected state about the axis 
£ by an angle given by Xxf 

exp{iA%^'}|Tx.), = (29) 

Thus, measurement of the relative positions allows us to definitely determine how 
much expjiAxjiXe} rotates \\x)j (i-e. to Therefore, eq. [29]acts as a correction to 

the ensemble: instead of the original ensemble of pre-selected \^\ x ) and post-selected 
states, we will have a new ensemble with shifted pre- or post-selections. 
How could the relative positions be measured? Procedurally, we first measure 
the momentum of the photographic plate after the N particles have deposited their 
momentum. When we subtract from this the initial momentum of the photographic 
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plate, then we can determine the shift in the sum of the momentum for the N particles 
as a result of the WM interaction in a new way. After the final measurement of 
Ancn we then measure the N individual positions (i.e. Q J md ) of each particle that 
is deposited onto the photographic. Now, measurement of Pj^J will disturb the 
individual Q J md s but will not disturb the relative positions (since they commute with 
the total momenta). Therefore, even though the subsequent measurement of the TV 
Qmd s wm be un-related to the value of the Q^s during the WV, we can deduce what 
the relative positions were at the time of the WM through the individual positions. 

After substituting the single particle result eq. [2H] for the jth particle (i.e. using 
the rotated bra or ket), back into the N particle eq. [28], we have: 

II<Tyliexp{i-=g22^}l*>iexp - ^ (30) 

j=l ViV I 4 l A VmdJ J 

It is clear from eq. [30] that when we look at the particles that are left unknown 
after using all the information (both relative positions and the total momenta) and 
consider them as the final total spin, then it is like a robust experiment but now 
the coupling to each spin is and thus the remaining effect of the coupling in the 
exponential will be small. 



3.2 Proving the legitimacy of WVs for a new regime 

After using these corrections, we can now prove the validity of the WV approxima- 
tion. We will show that the final state of MD, i.e. of eq. [30] will be: 

$* n d = exp{^i V ~ai} exp{ } (31) 



(o"4 is the WV for the jth particle - a tilde will always refer to WVs calculated with 
rotated states) When this is transformed back to the momentum representation (as 
was done in the WV approximation used in eq. [TO]) , then the momentum of MD is 
shifted by the WV; i.e. the change in (the change in the sum of momentum 

N pi \ N 

SP^ = 6^%=-^^ai (32) 



i=1 VN VN i=1 

To prove the legitimacy of this WV calculation, we first assume for simplicity a small 
variance in the rotations so that each particle yields approximately the same WV, 
enabling us to re-write eq. EJJ as: 



\-n(N) \/n( N ) (0 md } 

, A Vmd , •- ■ A WmdlJV„ 4(Aq( n » 



(N)^2 

2 



{cos-^^ + ia w sin-^} N e 4(AQ n^ (33) 



5 If the uncertainty in the individual Q m d's is AQmd ~ lj then the uncertainty in Q^) is also 
AQ m d ~ 1 (because the spread in YlQmd is V^V and thus ^-^° d as 1). Therefore, the spread in 
^-^ md is negligible and thus Xi = Q l md — a ^ so nas tne same uncertainty as AQ m d ■ 
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Now, in order to perform a valid WV calculation, this function needs to be peaked 
around = 0. As long as there are no regions in which the size of eq. [33] (i-e. 
eq. exceed the exponential of MD then it will be as if we are around Qwl = 0- 
In other words, the legitimacy of the WV calculation can now be understood as a 
competition between the (scalar product) A term and the (probability) B term: 

1*1 = {l + ^-D^^f exp {"igfiy < 34 > 

A B 

If the quantity eq. [M] goes to for large Q then the WV approximation is valid. On 
the other hand, if the increase in A was not counter-balanced by the decline in B 
then we could not restrict the WV approximation around Q = because it would 
be much more likely to be located around large Q. Thus, the meaning of the new 
WV approximation presented here is that there is no other region in which the size 
of eq. [331 i- e - ec L- EH> is significant, except around Q = 0. We now ask what is the 
maximum value of A such that we still obtain a shift in the pointer by J2j=i <^w? 
We will see that the constraint A 1 as was required in §2.11 (and in the first article 
on WMs [2]) is unneccessary in order to obtain a valid WM. In fact, with RWM the 
coupling to each individual spin just needs to be A ~ 1. 



3.2.1 oo, Q finite: 



First we consider how large A can be for a legitimate WV in the regime N — ► oo 
with Q finite. In this case, we consider again the magnitude (eq. |34|) . A can be 

written as {1 + (a 2 — 1) }^~. If Q is finite when N — > oo then we can 

neglect and higher terms from the expansion of sin 2 . Thus A w {1 + (a^ — 
1 \^Ql^ ~ eX p| ( Q w-i) A2 Q 2 a s i ong as (q,2 _ -q^2 < ^ t j ien fty e increase in A, i.e. 



N 



exp{ ' Sjli is counter-balanced by the decline in the Gaussian B, and thus 
eq. [Mfis centered around Q = 0. 



3.2.2 Finite N 

We consider finite iV where there is a proper limit in which iV increases and the 
interaction goes to but A 3 is negligible. If we fix A and choose an iV such that 
XyN > 1, then we can measure the average exactly. The uncertainty of P m( j for 
N particles is \fN and the momentum grows as XNa > VN which implies that 
2\fNa > 1. Nevertheless iVA 3 is still small (i.e. NX > \/N) but iVA 3 < ^ and 
iVA 3 is the extra correction. For each spin there is a correction proportional to A 3 
which for iV particles is iVA 3 which is small compared to yN so \^NX 3 < can be 
neglected. We plot N = 20 (fig. HJ) to show that eq. [31] is an accurate approximation 
to eq. EDI 



18 








(d) 




Figure 4: Numerical comparison of the RWM wavefunction (a-b) with the ideal shift 
by the WV (045°)™ = 1/2 (c-d) for N = 20. a) Real part of eq. EOJ b) Imaginary 
part of eq. [30], c) Real part of eq. EH d) Imaginary part of eq. [3TJ 



3.2.3 MD shifted by WV 

Now X)j=i <^w ~ Na w and thus all N particles will deliver a momentum shift to 
the photographic plate of XNa w (and a shift in Pj^ of Xy^Na w ). The shift goes 
up as Na w while the uncertainty goes up as \/N (the variance is ((Apj^) 2 ) = 

[(AP md ) 2 + ((A(cr^)„) 2 )]) and thus ™ ( d N { « v^V. By choosing a sufficiently large 

> md / 

N, the single trial WM outcome can be arbitrarily amplified. We have thus shown 
that it is valid to perform a WV approximation (i.e. to replace eq. [30] with eq. |3T|) 
in a significantly stronger coupling regime, i.e. for A ~ 1, by measuring a variable 
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where it's shift is large compared to it's noise and thus, this is a RWM. 



3.2.4 Quantum Average of WVs 



In the last 2 regimes, we chose, for simplicity, to ignore the details of a significant 
variation in the WVs, e.g. 5"i ~ Na w . Even if the variance were significant, 

it is easy to see that our result is still valid. I.e. even if the individual \I/'s in 
the composite state are very different, then eq. [15] is still valid, only the average 
will be over different pre- and post-selections. While it is appropriate to replace a 
by it's WV, there are 2 reasons that rotations can be induced in the pre- and/or 
post-selection. Up to this point we have only discussed the first rotation (e.g. from 
\1x)j to which was corrected by the measurement of the relative positions. 

This produced a shift in the pointer by an average over well-known WVs. However, 

exp{ 1 ^ d Z)j=i can cause a second rotation in the pre- or post-selection if AAQmd 
given by eq. [33] is big enough. However, in the instant case, the shift generated in 
the pointer exp{i\Q^y/Na w } can be large, so even a small rotation to the pre- 
or post-selection will make a significant difference in the pointer shift. We leave 
this analysis of additional rotation in the pre- or post-selections which cannot be 
determined by measurement of the relative positions to a future article. However, 
we show here how the total momentum, , is shifted by a quantum average of 
WVs [7] due to the extra rotations by the additional uncertainty in Q ( after the 
relative-position corrections are made) with weights determined by the probability 
to obtain a particular Q that is associated with a particular WV as suggested by eq. 
[HI As a simple example, one may categorize the different WVs into different pre- or 
post-selections. Suppose a subset, rii, out of the ensemble of N particles will all be 
rotated to the same state (e.g. to |^i)) and thus will give one WV fj^, other subsets 
will be rotated to another state (e.g. to l^)) giving another WV fj^, etc. Using 
Y, n i = eq. [30] is re- written: 



exp{iAQ2 -7=774 } ' ' ' ex P{ iA< 5md- 



(N) n k _ 



fN 



exp < — 



4(AQ2) 2 



(35) 



3.2.5 Obtaining EWVs instead of just an ordinary WV 

. In order to implement the amplification scheme, we need to obtain EWVs. To 
obtain an EWV (i.e. outside the eigenvalue spectrum), rather than an ordinary WV, 
we need to control the rotation of the pre- or post-selection by controlling XAQ. E.g. 
if Q is limited (e.g. AQmd ~ 1) and A is limited to a particular range sufficient to 
deliver a EWV at every point of Q within AQ m d, then we will also obviously obtain a 
EWV for the quantum average of WVs and do not need to be concerned with other 
issues such as the slope of Q. However, anytime there is a way to get inside the 
spectrum of eigenvalues, there will be an exponential increase in the probability to 
obtain that WV. This can be seen from eq. M in that the fluctuation in the system is 
also relevant for the probability to obtain different post-selections: as the fluctuation 

6 For any product state, we still have that eq. [5] is exactly true but with an a = ^f?„, i.e. 

C^ (N) i* (N) > = Cd ) E^i^ N )>+Cd ) ^i* (N) ■!■>• 
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in the system increases, the probability of a rare or eccentric post-selection also 
increases. However an attempt to see this through WMs will require the spread in 
the MD to be increased and this increases the probability of seeing the strange result 
as an error of the MD. 

4 Conclusion 

In this article, we have introduced a new WM procedure for finite samples which 
yields accurate WVs that are outside the range of eigenvalues and which does not 
require an exponentially rare ensemble. This procedure was motivated by an ap- 
plication to quantum metrology which provides a unique advantage over the usual 
(pre-selected-only) approach in the amplification of small non-random signals: if the 
coupling A between system and MD is unknown and contains additional small er- 
rors that are not random, then actually performing a WM which yields an EWV 
can provide new information, e.g. by allowing us to distinguish the shift from the 
non-random force (incorporated into A) from the large EWV shift due to the WM 
interaction. The usual WV approximation becomes more and more precise in the 
idealized weak limit of AAQmd — > and N — > oo in which there is no disturbance or 
back-reaction on the system. Neither of these limits are realistic in practical ampli- 
fications because first of all, we must have a finite N, and second of all with a finite 
N, we must increase A to discern the WV from the noise. RWM minimizes the two 
uncertainties in determining the WV which arise due to 1) the inability of MD to 
measure definite WVs due to the MD's uncertainty AP m( j and 2) the back-reaction 
on the system due to AQ m d creates an uncertainty in the WV of the system itself 
By providing additional corrections for these uncertainties, RWM extends the cou- 
pling constant regime which thereby extends the potential utility of amplification of 
unknown forces [TTJ [8] . 

The RWM can also be used to augment the SWM of §2.11 given a large ensemble 
(N — > oo) because we can now interpret what is the average of WVs corresponding 
to the total momentum shift for a stronger coupling constant by calculating the 
distribution of pre- and post-selections through the distribution of relative Q's. We 
can calculate what the distribution in Q will be for N — > oo, even if we do not know 
the distribution for the individual Q's for any individual particles (the distribution 
of J2 Q becomes a Gaussian for large N for almost any individual distribution of Q) . 
However, [8] for finite N we cannot simply use a calculation because the fluctuation of 
the relative Q's becomes important and can only be obtained through measurement. 

Normally a valid WV calculation requires MD to be centered around = 0. 
However, in [7] ideal measurements were converted to WMs by post-selecting MD to 
be in a certain region of Q and in P m d, i.e. different regions of Q were sampled by 
multiplying by a function centered at Q CO m- i-e. a function of Q' = Q — Q CO m such 
as exp ~^ Q r^2° n ^ which is like starting the MD not with Q = but with Q = Q com . 
Results centered at different Q = Q CO m are then summed. However, even such limited 
projections can still disturb each other. The new RWM presented here is more subtle 
because the relative coordinates commute with the total momentum and so can be 
simultaneously measured without disturbing each other. By measuring the relative 
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positions, we can go beyond the weak approximations used in the past (i.e. A< 1). 
Since we are able to measure the relative positions exactly, we are also able to make 
these corrections exactly. We thus have a much stronger interaction (i.e. a A that 
does not have to be <C 1) and still we can obtain EWVs. 

Acknowledgments: The author thanks Yakir Aharonov and Alonso Botero for 
many fascinating discussions. 
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